Origin of spontaneous electric dipoles in homonuclear niobium clusters 
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Surprisingly large spontaneous electric dipole moments recently observed in homonuclear niobium 
clusters below 100 K (Moro et al. 300, 1265 (2003)) are explained using first-principles electronic 
structure calculations. The calculated moments for Nb n (n < 15) closely follow the experimental 
data in which large dipole moments are seen for n = 11-14. We establish that the dipoles are 
strongly correlated with the geometrical asymmetry of the clusters. The magnitude of the dipole 
moment is roughly proportional to the spread in the principal moments of inertia and its direction 
tends to align along the axis of the largest principal moment. Charge deformation densities reveal 
directional, partially covalent bonds that enhance the formation of asymmetric geometries. Classical 
simulations of the deflection of a cluster in a molecular beam reveal that the electronic dipole may 
persist at higher temperatures, but is masked by the rotational dynamics of the cluster. 



Though the response of nanoparticles to an external 
electric field has been actively studied for several years, 
the recent observation of permanent electric dipoles on 
the order of several debye (D) in free Nb„ clusters 
(n ~10-100) has introduced an entirely new and unex- 
pected dimension. That the intrinsic moment is only 
observed at temperatures T(n ~ 10) « 100 K[l| has led 
to speculation that the origin of large electric dipoles 
might be related to superconductivity in bulk Nb (rather 
robust at T c =9 K) and that the dependence on temper- 
ature may signal a nanoscale "ferroelectric" transition. 
Permanent dipoles in mixed clusters and molecules, in- 
cluding canonical polar molecules like H2O and HF, are 
well understood in terms of charge transfer between un- 
like atoms. However, it is more challenging to envision 
what forces could drive the underlying charge separation 
in homonuclear clusters of metal atoms. This unexpected 
and nonintuitive dielectric behavior in Nb clusters [l| has 
also attracted attention to the conceptually related ques- 
tion of quantum dipoles in symmetric clusters[ll3. 

A fundamental requirement for the occurrence of a per- 
manent dipole is the lack of a center of inversion. For a 
homonuclear system, however, the issue is more complex 
because disproportionation into charged ions is energeti- 
cally disfavored. These questions are more perplexing in 
clusters that are going to become metals (which support 
no charge separation) as their size increases. The calcu- 
lations we discuss demonstrate that the origin of these 
large dipoles does not require a new state of matter re- 
lated to superconductivity as suggested in Ref. fl, but is 
explained directly in terms of directional partially cova- 
lent bonds that favor the formation of clusters with low 
symmetry. The resulting charge separation, which im- 
plies internal electric fields on the order of 10 6 V/cmQ, 
is sustained by the chemical forces of the Nb bonding. 
By performing classical simulations of the rotational dy- 
namics of the clusters in a molecular beam, we provide 



evidence that thermal averaging, instead of a loss of in- 
trinsic moment, accounts for the observed temperature 
dependence. 

The atomic and electronic structures of neutral Nb n 
clusters (n < 23) have been extensively studied using 
first- principles electronic structure met hods 0, IH 
but there has been no theoretical study of their dielec- 
tric properties. In this Letter, the electric dipole mo- 
ments of Nb„ clusters (n < 15) are analyzed in terms 
of their asymmetry and charge deformation using the 
electronic structure codes Gaussian980 and AbinitjjJ. 
Calculations with Gaussian98 used Stuttgart-Dresden ef- 
fective core potentials^ and the B3PW9l[lll E3 hy- 
brid exchange-correlation functional. Calculations with 
Abinit0 used the planewave (norm-conserving) pseu- 
dopotential method 20] and PBE960 GGA functional. 
The lowest energy structures as well as isomers close 
in energy (within 25 meV/atom) obtained by Kumar 
and Kawazoe0 were used, except for Nbi3, for which 
a new structure lower in energy was found. Isomers of 
Nb+ clusters have been observed in spectroscopic and 
reactivitv|l5| experiments. 

The calculated permanent electric dipoles (TABLE P) 
are compared in FIG. ^with the reported valuesQ for 
Nb clusters emitted after reaching thermal equilibrium 
with He gas at 50 K. Considering the many uncertainties 
in presenting such a comparison, the level of agreement is 
satisfying: moderate moments are found for 3 < n < 9, 
large moments for 11 < n < 14, and essentially zero 
moment for n = 4, 10, and 15. The only notable dis- 
crepancy is Nbin for which experiments show significant 
dipole moment, |2lj while theoretically it is zero because 
of the symmetric nature of the lowest energy structure. A 
general feature in FIG. 1 is the larger calculated dipole 
moment in comparison to experiment. This difference 
may be due to the neglect of rotational dynamics in the 
experimental analysis[lj. Calculations on isomers, such 



1 1 Exp 
la -A Theory 




FIG. 1: Comparison of the calculated dipole moments and 
"low temperature" (50 K) experimental data from FIG. 2 of 
Ref. (1 D = 0.20819 eA). Higher energy structures are 
labeled in alphabetical order, following the notation of Ref. 
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TABLE I: Principal moments of inertia I1-I3 (MNbA 2 ) and 
dipole moment /j, (D) for Nb n - 
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as Nbt)& and Nbub, whose binding energies are lower by 
1 and 9 meV/atom respectively, show dipole moments 
closer to the experimental data[l|, which may be an in- 
dication of isomers in the experimental conditions. How- 
ever, Nb6a is in good agreement with experiment in com- 
parison to Nbgfc, whose binding energy is also only mod- 
erately lower (20 meV/atom). 

For n > 2 only in the case of Nb4, which is a regular 
tetrahedron^, does symmetry strictly forbid an electric 
dipole. Other cases of vanishing dipole moment, such as 
Nbio and Nbi5, are attributed to near symmetries. Nbio 
has an approximate symmetry, whereas Nbis has 
near inversion symmetry with a cubic structure. These 
symmetries are reflected in the cluster's principal mo- 
ments of inertia 7i, h, and ^3 (TABLE P). Nbzi, for 
example, is a spherical top {I\ = I2 = I3), and Nbio and 
Nbis are symmetrical tops (I2 — I3 / h), although for 
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FIG. 2: The strong correlation between the dipole moment 
(green, right axis) and asymmetry of the cluster, as quantified 
by the difference between the largest and smallest principal 
moments of inertia (blue, left axis). 



Nbis the deviation from a spherical top is quite small. 
The rest of the clusters studied are asymmetrical tops 
(h 7^ I2 7^ h)- One measure of the degree of asymmetry 
is the spread A J = J max - ^ mm) where I max (J mm ) is 
the largest (smallest) principal moment of inertia. Asym- 
metric clusters have correspondingly large values of AI. 
FIG. |2] shows the strong correlation between AI and the 
calculated electric dipole. Since the principal moments 
of inertia map the dynamics of a rigid body to that of an 
ellipsoid, AI can be interpreted as quantifying the exis- 
tence of a preferred structural axis in the cluster. This 
interpretation is supported by the alignment of the elec- 
tric dipole. In cases where AI is large (Nbii-Nbi.4), the 
direction of the dipole moment and the principal axis 
corresponding to /max are close to collinear. The 26° 
deviation from collinearity of Nbi2 is the only exception 
to this otherwise near-perfect (within 5°) collinearity. 

To understand the origin of the electric dipole we con- 
sider the charge deformation density due to bonding, that 
is, the difference between the charge density of the clus- 
ter and that of the isostructural spherical neutral atoms, 
for which the dipole moment vanishes. FIG. E3 shows the 
charge deformation for the simple and illustrative case of 
Nb^|22j. an isosceles triangle|j| with a dipole moment of 
0.5 D. In FIG. |3 two isosurfaces are shown. The blue 
(red) surface corresponds to a negative (positive) isocon- 
tour of value -2.0 e/A 3 (+2.0 e/A 3 ). Taken together, 
these isosurfaces show how charge redistributes during 
the formation of bonds. Two general features can be dis- 
cerned. First is the formation of a few specific covalent 
bonds as revealed by the charge build up between ions, 
which reflects charge being pulled in from elsewhere in 
the cluster. Second, some charge is pushed outward at 
each surface site due to Coulomb repulsion. 

FIG. 0] shows the charge deformation for Nbi2 and 
Nbi5, which are representative of clusters with large 
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FIG. 3: Charge deformation of ND3 (0.5 D). Due to the mirror 
symmetry, the dipole moment (orange arrow) is restricted to 
lie on the symmetry axis. 



dipole moment and no dipole moment respectively. The 
same color scheme and isocontour values are used as in 
FIG. |21 and the general features remain. Only in cases 
with a dipole moment, such as Nb 12 , is there a strong 
asymmetric character to the deformation density. This 
reveals how strong directional bonding in isolated regions 
of the cluster leads to the intrinsic charge separation 
giving rise to asymmetry and the permanent dipole mo- 
ment. We find however that the regions in the cluster in 
which the charge deformation shows this additional co- 
valent character don't always coincide with the shortest 
bonds. Instead, directional bonding between rf-electrons 
is important, driving the development of low-symmetry 
geometries. 

An additional aspect of the experiments of Moro et 
al. is the disappearance of the dipole (i.e., a change in 
the deflection of the cluster beam) when the clusters are 
formed at higher temperatures 1], which has been likened 
to a ferroelectric transition. One possible mechanism for 
the disappearance of the polarized electric state is the 
thermal excitation of vibrational modes within the clus- 
ter, which should be directly related to the temperature 
at which the clusters are formed. The charge response 
to vibrations might disturb the intrinsic charge separa- 
tion in the static cluster at equilibrium. In addition, if 
the direction of the dipole moment were to fluctuate, this 
charge separation, if still in existence, could become un- 
observable if the time average of the dipole moment (/2) 
were sufficiently reduced. Such vibration-induced fluctu- 
ations of the dipole moment are known to be important 
in weakly bound molecules 0. We suspect that such 
vibrations are responsible for the zero dipole moment 
observed in Nbs and Nb7 which have distorted trigonal 
and pentagonal bipyramid structures respectively. In the 
symmetric structure of M>7, the highest occupied molec- 
ular orbital is doubly degenerate but occupied with one 
electron, which leads to a Jahn- Teller distortion. When 




FIG. 4: Charge deformation for (a) Nb i2 (2.3 D) and (b) 
Nbis (no dipole). The direction of the dipole moment is rep- 
resented by the orange arrow. Note the formation of strong 
covalent bonding. Unlike systems with a large dipole moment, 
no strong asymmetric covalency is seen in Nbis. Instead the 
charge deformation density is entirely intra-atomic. 

formed and selected at a higher temperature, thermal 
motion may move the cluster through several structures 
with differing dipoles, the net effect of which would pro- 
duce canceling deflections. 

To quantify the influence of vibrations, we calculated 
the vibrational spectrum for Nbi2, a cluster with a large 
permanent electric dipole (2.3 D), and looked at the exci- 
tation of the 6 lowest energy vibrational modes between 
60-105 cm^ 1 . These are in the range of those that would 
be excited during the experimentally observed transition, 
which occurs at a temperature on the order of 100 K. 
The atomic displacements associated with these vibra- 
tions are less than 0.1 A and are unable to significantly 
alter the charge density or the permanent dipole, with 
relative changes of less than 3% or 0.1 D in magnitude 
and a few degrees in direction. This implies that an- 
other mechanism, beyond harmonic vibrations, is needed 
to explain the disappearance of the dipole moment. 

Classical simulations by Dugourd et aZ.^J, using the 
thermal occupation of rotational states, were able to re- 
produce the molecular beam deflection of TiCeo clusters, 
but full temperature dependence was not studied. FIG. 
shows the classical deflection of a representative nio- 
bium cluster (Nbi2) with a dipole moment of 2 D at 20 
and 300 K when the homogenous electric field is on and 
off. These profiles were obtained by numerically evaluat- 
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FIG. 5: Classical deflection of a model Nbi2 cluster with a 
dipole moment (i = 2 D at 300 K (top) and 20 K (bottom). 
At higher temperatures, the rotational dynamics masks the 
permanent dipole moment. 



effects in experiment accurately. However, because the 
structure of the deflection profile is richer at low tem- 
peratures, and depends directly on the inertial moments 
of the cluster, future low-temperature molecular beam 
experiments have the potential to measure structural pa- 
rameters that are directly comparable to theoretical pre- 
dictions. 
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ing the deflection of 10 6 random initial configurations 
(orientation and rotational energy) of a cluster with sym- 
metric inertial moments Ii = I2 = 45 M^ A 2 and I3 = 
65 MNb A 2 . The deflection due to the induced dipole, 
which would shift each profile a constant distance, has 
been ignored by setting the electric polarization a = 0. 
In addition, the broadening of the profile due to the clus- 
ter's reduced translational velocity at lower temperature 
has been divided out. 

The classical deflection profile FIG. |5] and the experi- 
mental profiles in Ref. Q] agree qualitatively. In particu- 
lar, the asymmetric shape at lower temperatures, which is 
an experimental signature of the permanent moment 0, 
is clearly seen in the profile at 20 K. At 300 K the de- 
flection profile is more symmetric, which has been inter- 
preted as a loss of the intrinsic moment but is actually 
the effect of additional thermal averaging. 

Finally we return to the question of (a lack of) sym- 
metry that we began with. Our calculations have traced 
the origin of large spontaneous electric dipoles to the 
structural asymmetry of Nb„ clusters that is enhanced 
by the directional bonding favored by d electrons. These 
correlations can be sharpened. First, the large asym- 
metric inertial moments observed in systems with large 
electric dipole moment result from oblate, not prolate, 
cluster shapes. Second, although the magnitude and 
direction of the electric dipole correlates very strongly 
with an (asymmetric) imbalance of inertial moments, 
this cannot be the whole story, since breaking reflection 
symmetry — of which the inertial moment themselves pro- 
vide no information — is required for dipole formation. 

That the apparent temperature dependence of the elec- 
tric dipole can be explained classically implies a more 
rigorous treatment of rotational dynamics 18] may be re- 
quired to distinguish induced versus permanent dipole 
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